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We consider Mach-Zehnder and Hong-Ou-Mandel interferometers with nonclassical states of light 
as input, and study the effect that dispersion inside the interferometer has on the sensitivity of phase 
measurements. We study in detail a number of different one- and two-photon input states, including 
Fock, dual Fock, NOON states, and photon pairs from parametric downconversion. Assuming there 
is a phase shift <?!>o in one arm of the interferometer, we compute the probabilities of measurement 
outcomes as a function of 4>o, and then compute the Shannon mutual information between 00 and 
the measurements. This provides a means of quantitatively comparing the utility of various input 
states for determining the phase in the presence of dispersion. In addition, we consider a simplified 
model of parametric downconversion for which probabilities can be explicitly computed analytically, 
and which serves as a limiting case of the more realistic downconversion model. 



PACS numbers: 42.50.St,42.50.Dv,07.60.Ly,03.67.-a 

I. INTRODUCTION 

Interferometry is both an important tool for practical 
measurements and a useful testing ground for fundamen- 
tal physical principles. As a result, the search for meth- 
ods to improve the resolution of interferometers forms an 
active area of study. It has been shown by a number of 
authors ([E [3, S, 01) ^^^.t nonclassical states, in particu- 
lar those with high degrees of entanglement, when used 
as input to an interferometer can lead to resolutions that 
approach the Heisenberg limit, the fundamental physical 
limit imposed by the uncertainty principle. Most of this 
previous work has dealt with idealized interferometers, 
with no dispersion or photon losses. Before quantum 
interferometry may become a useful practical tool the 
question must be asked as to how well the conclusions of 
these previous studies hold up in more realistic and less 
idealized situations. In this paper, we will attempt to 
take the next step along this road by adding dispersion 
to the apparatus and examining what effect this has on 
the phase sensitivity of interferometry with nonclassical 
input. The motivation for this work is the desire to ulti- 
mately construct quantum sensors that can measure the 
values of external fields by measuring the phases shifts 
they produce in an interferometer. 

In particular, the nonclassical input states we will con- 
sider are (i) Fock states | A'', 0) which have a fixed number 
of photons incident on one input port, (ii) dual or twin 
Fock states | A^, A^) which have the same number of pho- 
tons incident on each input port, and (iii) NOON states 
^ [\N, 0) + |0, N)]. Here, |A"a, Nb) denotes a state with 
numbers Na and Ni, of photons entering each of the two 
interferometer input ports. 

There has been a great deal of recent work on the pro- 
duction of nonclassical states of light with large (A^ > 2) 
numbers of photons by means of postselection (for exam- 



ple, [1, [1, 0, Q), however at present the utility of these 
postselection schemes for application to practical situa- 
tions is not clear. Although this work is useful for clar- 
ifying the scientific issues involved, it is not technologi- 
cally feasible at present to use these methods to produce 
the desired states on demand. Rather, postselection pro- 
duces states statistically, at random times, and therefore 
can not be relied upon to produce states on demand for 
a quantum sensor. In addition, for large photon number, 
great care must be taken to distinguish between states of 
A^ photons and those of A^ — 1 photons, making it diffi- 
cult to prevent mixed states from appearing, which would 
change the physics involved. In contrast, two-photon en- 
tangled states with well-defined properties can be easily 
produced by parametric downconversion or other meth- 
ods. 

Due to the current practical difficulties of producing on 
demand entangled photon states with large, well-defined 
N , we save the large A^ case for later study and restrict 
ourselves in this current paper to situations which are 
both simpler and of more immediate practical interest, 
namely the cases of one or two photons. Furthermore, 
for the two-photon case, we consider two possibilities: (i) 
the photons may be uncorrelated in frequency, or (ii) the 
pair may be produced through spontaneous parametric 
downconversion (SPDC), resulting in anticorrelation be- 
tween the two frequencies. 

Our goal is to compare the usefulness of each of these 
cases for making phase measurements in the presence of 
dispersion, so we will need a means of quantifying the sen- 
sitivity of the interferometer with respect to these mea- 
surements. Consider a single shot consisting of a non- 
classical state of light with a fixed number of photons 
being injected into the input ports of the interferometer. 
Suppose some phase-dependent observable M[(j)) is mea- 
sured during this shot. The usual way to define the phase 
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sensitivity of the measurement is by computing 

-1 



A0 



AM. 



(1) 



However, this is correct only if the probabihty distribu- 
tion of the phases has a single peak and is approximately 
Gaussian in shape. A more general strategy is to take an 
information-theoretical approach and to define the quan- 
tum fidelity by means of the Shannon mutual information 
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Here, m and are the measured values of the random 
variables M and $, while P{m\(t>) is the conditional prob- 
ability of obtaining measurement m given the phase (p 
on a particular shot. In this formula, we have also as- 
sumed maximum ignorance of the phase, i.e., we have 
assumed a uniform distribution for 0, p{(j)) = Sup- 
pose that the detectors have a characteristic time-scale 
Tq. Then in this context, a single shot will consist of 
a well-defined number of photons entering the appara- 
tus simultaneously (i.e. within a temporal window much 
smaller than T^) and seperated in time from any other 
entering photons by a time > Td- The mutual infor- 
mation is a measure of the information gained per shot 
about the phase $ from a measurement of the observ- 
able M . In our case, the role of M will be played by 
the number of photons detected at each of the output 
ports. For N input photons, output detector C will count 
I photons, detector D will detect the remaining N — I 
photons, and the sum in equation [2] will become a sum 
over where / — 0, 1, . . . iV. Throughout this paper we 
will use the quantum fidelity as our measure of phase 
sensitivity. Besides being of very general applicability 
and giving a precise, calculable measure for the utility of 
a measurement, the introduction of the mutual informa- 
tion provides a link to the theory of quantum information 
processing. Bahder and Lopata have computed the 
quantum fidelity as a function of N for idealized lossless 
and dispersionless interferometers with Fock and NOON 
state inputs. In the following sections, we will see how 
their results change for the cases of iV = 1 and N — 2 
when dispersion is present. 

Although not the principal focus of this paper, it 
should be noted that the existence of multiple peaks 
in the output probability distributions invalidate the as- 
sumptions used to derive the Heisenberg bound from the 
Cramer-Rao lower bound, which makes input states with 
multimodal distributions especially interesting from the 
point of view of the study of phase sensitivity. Note 
that violations of the Heisenberg limit have recently been 
shown to exist in another context, distinct from the sit- 
uation examined in this paper, namely in the context of 
nonlinear interferometry [HI [HI . 

We will assume one branch of the interferometer has a 
dispersive element which gives the photon wavenumber 



k a frequency dependence of the form 

k{uj) — ko + a{u} ~ bJo) + f3{i^ — tjJo)'^ , (3) 

ignoring the possibility of higher order terms. The other 
interferometer arm will be assumed to be of negligible 
dispersion. Here, a is the inverse of the group velocity, 
and P is the group delay dispersion per unit length. 

In addition to the Mach-Zehnder interferometer, we 
will examine the fidelity of an alternate setup used in 
[Q], in which NOON states are incident on a single beam- 
splitter used as a Hong-Ou-Mandel (HOM) interferome- 
ter. We will then be in a position to compare the pos- 
sible input states and interferometer setups, with a view 
to gaining insight into their relative usefulness in prac- 
tical measurements. In the two-photon cases, we must 
distinguish between situations in which the photon en- 
ergies (or frequencies) are correlated and those in which 
they are independent. Thus, after we examine the case of 
energy-uncorrelated photons, we look at photon pairs an- 
ticorrelated in energy. We further consider two subcases 
of the latter: (i) a simple model which can be solved 
analytically and which amounts to a simplified version 
of spontaneous downconversion, and (ii) a more realistic 
but less analytically tractable version of downconversion. 

The plan of this paper is as follows: in section |TT] we 
consider the setup for the dispersive Mach-Zehnder inter- 
ferometer and define the input states we will use in more 
detail. In section Hill we apply the possible one-photon 
inputs to the interferometer and compute the probabili- 
ties for the various possible outcomes. In sections HVl and 
IVl respectively, we do the same for the Mach-Zehnder 
interferometer with several different two-photon inputs 
and for the HOM interferometer with N = 2 NOON state 
input. In section IVII we compute and plot the mutual 
information for each of the preceeding cases as functions 
of bandwidth and dispersion levels; we then compare and 
discuss the results for the various cases. Finally, in sec- 
tion lVIll we repeat the same calculation for input consist- 
ing of a photon pair produced via spontaneous paramet- 
ric downconversion before arriving at final conclusions in 
section IVIIII 

For ease of reference later, table I summarizes the spe- 
cific cases we will examine over the following sections. 



II. THE DISPERSIVE MACH-ZEHNDER 
INTERFEROMETER 

Consider the Mach-Zehnder interferometer of figure (TJ 
with 50/50 beamsplitters. Assume for the moment that 
there is no dispersion in the apparatus. Let a^j and b^^ be 
operators that annihilate photon states in the two input 
ports A and B. They obey the usual canonical commuta- 
tion relations with the corresponding creation operators 
al and bj- 
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TABLE I: Summary of the special cases examined in the later 
sections of this paper. 



Case 


#of 


Interferometer 


Input 


Frequency 


No. 


photons 


Type 


State 


Correlation 


A 


1 


MZ 


Fock 


not applicable 


B 


1 


MZ 


NOON 


not applicable 


c 


2 


MZ 


Fock 


none 


D 


2 


MZ 


Dual Fock 


none 


E 


2 


MZ 


NOON 


none 


F 


2 


MZ 


Fock 


anticorrelated 


G 


2 


MZ 


NOON 


anticorrelated 


H 


2 


MZ 


Dual Fock 


anticorrelated 


I 


1 


HOM 


NOON 


none 


J 


2 


HOM 


NOON 


none 


K 


2 


HOM 


NOON 


anticorrelated 


L 


2 


MZ 


SPDC Fock 


anticorrelated 




FIG. 1: Mach-Zehnder interferometer with dispersion in one 
arm. There is also a phase shift (po of nondispersive origin in 
the same arm. 



with all other commutators vanishing. For independent 
photons, the input states to the interferometer can be 
described in terms of the number of photons entering the 
two ports: 



\Na,Nb;uJi,. 



At 



a. 



|0), (5) 



where Na and Nb are the number of photons in ports 
A and B, respectively, and |0) is the vacuum state with 
no photons. Similarly, Nc, Nd, c^^, and d^^ will represent 
the photon numbers and annihilation operators at output 
ports C and D. 

The effect of the Mach-Zehnder interferometer on a 
given input state may be described in terms of the scat- 
tering matrix, S{(j)). The initial and final annihilation 



operators are related by a scattering matrix S{4'): 



(6) 



where is the relative phase difference experienced by 
photons in the two arms. In the absence of photon 
losses in the system, the scattering matrix will be uni- 
tary. Then, for an ideal Mach-Zehnder interferometer, 
the scattering matrix is given by 



5(0) 



where the Pauli matrices are 



and 




1 
-1 



(7) 



(8) 



In this scattering matrix we have assumed (as we will 
assume henceforth) that the lengths of the two interfer- 
ometer arms are equal, Li = L2- Using this matrix in 
equation [SI we can invert the equation and take adjoints 
to arrive at the following result: 



a!, sin 



— — a ' cos — 
2 2 



cl cos -j^+dl sm - 



.i4>/2 



(9) 
(10) 



We assume that the frequency distribution for each 
incoming photon is Gaussian and that each Gaussian 
has the same width and central frequency, of the form 

1 / ■,2 

g-^cr(ui-uio) _ Input and output states will either be 
states of definite photon number in the sense that they 
are eigenstates of number operators of the form Nj = 

J (kua'^J^^ a^^ (where a^^ is the annihilation operator for 
photons at the jth port), or else superpositions of such 
states. 

We introduce dispersion to the upper branch of the 
interferometer by giving the wavenumber k a frequency 
dependence of the form in equation [3) We assume that 
the dispersion in the other branch of the interferometer 
is negligible, i.e. that k{oj) = kg in that branch. The 
length of the portion of the upper arm for which fc((jj) 
differs from fco will be denoted L, where < L < Li. In 
addition to any phase difference resulting from the asym- 
metric dispersion, we also assume that photons travel- 
ling through the upper branch of the interferometer gain 
an additional phase difference (po relative to the lower 
branch. (j)o is any phase difference of nondispersive ori- 
gin that may be present in the setup; this may be due to 
a difference in path length, or an interaction of one arm 
of the interferometer with an external field. Note that for 
our setup, the assumption of a balanced interferometer 
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entails no loss of generality; to account for an unbalanced 
interferometer, it suffices to simply include a term of the 
form ko{Li — L2) inside the phase factor 0o- 

In the presence of the dispersion, the scattering matrix 
will now be of the form 



S{4>o) 



— _jg«'=oiigj0(t-')/2 



e«0('^) _ 1 -i (e»"^'('^) + 1) 
i (e*'^('^) + 1) - (e''"^('^) - 1) 
0M " 



(11) 



where for future convenience we have shifted the fre- 
quency dependence into a new phase angle by defining 



{lo) = 00 + aL{uo - Wo) + l3L{uj — wq)^ 



(12) 



Consider N photons entering the interferometer and 
assume for now that their frequencies are independent 
variables. The Fock, dual Fock, and NOON input states 
are of the form: 

x4,---«L|0) (13) 
^i.jO> (14) 



X a], . . . aj, b}, 



and 



x/2 



[|7V,0) + |0,7V)] 



1 1 ^(J\N/i 



X 



|o), 



(15) 



where the bandwidth of the incident beams is given by 
Auj = a^^l"^ . If the photons are produced by SPDC, 



then the frequencies must occur in pairs with the pho- 
tons in each pair being equal distances above or below 
the pump frequency; we will consider this situation in 
simplified form in section IIVBI and in a more realistic 
form in section IVTIl 

Suppose that one of the A^-photon or 2A^-photon states 
described above is input to the interferometer. Write this 
input state as iV'in)- Then, assuming that the frequencies 
of the final photons are not measured, we want the joint 
probabilities to find iVc photons at detector C and Nd 
photons at detector D (with A^c + = Na + Ni,) for a 
given nondispersive phase shift (po in the upper interfer- 
ometer arm. These probabilities can be expressed in the 
form 

PiN,,NM = {^^n\^T{N,,Nd■,(|)omn), (16) 
where the projective operator Tt{Nc, Nd] 0o) is defined as 



7r(7V„ iVrf, 0o) = / dn\N,,Nd;n,(f>o){N,,Nd;n,(f>o\ 



with 



\N,,Nd;n,cl)o) 



(17) 



VNJ Nd\ 



i £ , 



JO). 
(18) 

Here we have suppressed the 0o-dependence of the 
Co and da' operators for notational simplicity, and 
have represented the collection of output frequencies 
{rii, . . . , . . . , by the single symbol il. Sim- 

ilarly, dfl is being used as shorthand for the full frequency 
integration measure dVli . . . dD,j\j^dCl[ . . . dQ'j^^. These 
probabilities may also be expressed in the form 



P{N,,Nd\<l>o)^ / dn\{N,,Nd■,n,<j>o\^p^n)f 



(19) 



From equation [2] the mutual information between the 
phase $ and the output photon numbers M is then 



: M) 



27r ^ 



d(t>o P{N,,Nd\M^og 



2TrPiN,, Ndl^o) 
/:^d0o P{N,,Nd\M 



(20) 



We note that the probabilities P{Nc, Nd,\4>o) are also 
conditional upon the values of a, /?, and a, although we 
do not explicitly include these parameters in the notation 
for the probabilities for the sake of notational simplicity. 
We now restrict ourselves to the cases iV = 1 and N — 
2, and proceed in the following sections to compute the 
mutual information H for a number of different possible 
input states. 



III. MZ INTERFEROMETRY WITH 
ONE-PHOTON INPUT 



In this section, we begin with the cases in which there 
is only one photon in the initial state. 
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Case A: One-photon Fock state. We introduce the 
normalized input state 



\^,^), ^ |10), - (/- / duj e-^^(— '')'ai|0), (21) 



representing a single photon incident on port A. Using 
relations [51 and [TUl this is equivalent to 



|10)c 



2 y TV 



[cl (e'^("' - l) - idl (e'"^'"' + l)] |0). (22) 




The photon may leave the interferometer via either 
port C or port D. We assume that the detectors count 
the number of photons leaving the apparatus but do not 
measure their frequencies. Therefore, we must integrate 
over the final frequencies. The output state is then mea- 
sured using the projective operators 



^(1,0) = Jdncl 



CO 



(23) 
(24) 



Expectation values of these operators give the proba- 
bilities of measurement outcomes: 



^(0,1) = / dndl\0){Q\dn. 



1 + 



^ COS N-o + -- 



(25) 



(26) 



In the previous two lines, we have introduced some 
notation that will be convenient for simplifying the re- 
sults of this and the following sections. The parameters 
ri,r2, 61,62 are defined by (see figure [2]): 



^ 1 + 1^^ 

cr 



(3Ly 

2a J 



tan 01 — 



tan 60 



a 

2^' 



(27) 
(28) 



Note that these parameters depend on the second order 
dispersion coefficient (3, but not on a, and that when (3 
vanishes we then have ri — r2 ~ \ and 61 = 62 = 0. 
Case B: One-photon NOON state. The input state 

is 



1 



where 



V2 



101)] 



J 



{cl[i 
+ dl 



duj e-t^^-'^o)' 



2V2 



(29) 



- 1) - (^ + l)e'''^(") (30) 



FIG. 2: Definitions of n, r2, 9i, and 82. 



The resulting output probabilities in this case turn out 
to be 



1 ^— sin 00 + — - ^ 2 2 

2 Arfa-' 



1 -I =— sin 00 + — - ^ 2 2 

2 Arfa-' 



(31) 



(.32) 



IV. 



MZ INTERFEROMETRY WITH 
TWO-PHOTON INPUT 



We now consider input states with two photons distributed in as- 
sorted ways among the input ports. However, now we must make 
a distinction as to whether the two photon frequencies are inde- 
pendent or correlated in some manner. We treat the uncorrelated 
version first. Then we will examine one particular case of frequency- 
correlated photons which is of special interest for experiment: that 
of photon pairs created through spontaneous parametric downcon- 
version (SPDC). In this section we treat only a simplified version 
of SPDC which will allow us to obtain simple exact expressions for 
the probabilities of all of the output states. In a later section we 
will compare this simplified SPDC to a more realistic version for 
which only numerical results are available. 



A. Two-Photon Input with Uncorrelated Energies 



Case C: Energy- uncorrelated two-photon Fock state. 

Sending a two-particle Fock state into input A, 



12, 0)<T = \/ ^ / duidu2 e 2 



[(cJl-t^0)^ + ("2-"0)^]st -t 



(33) 

where we use equations l9land ll0l to write each a) factor in terms of 
the output operators and . After a straightforward calculation, 
this leads to the following output probabilities: 
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P(2,0|</.o) 



1 cPpL^ 
ri V 2 4rf(72 



(34) 



P{0,2\4,o) 



f(l,l|</'o) 



1 H —— cos <po + — - 

'n V 2 



1 cos <P0 + , , , 



(35) 



(36) 



Case D: Energy-uncorrelated two-photon dual Fock in- 
put. The normalized input state is 



4 V T J 

X {-icl.ci^ (e'^(-i) - 1) (e"^("^) + l)+idi^di^ - l) (e'^(-i) + l) 

(e"*''*'^' - 1) (e'^(-2) _ 1^ _ £t ^rft^^ + 1^ (^e'^t"^) + l) } |0>, 



(37) 
(38) 



which gives the results 
P(2,0) = P(0,2) 
1 



Ml ' 
2 I n 



21-'^ CT 



2(^0 



2(^0 + 6'i 



2r2cr2 



2r2a2 



(39) 



(40) 



Case E: Energy-uncorrelated two-photon NOON state. 

For the input state 



spontaneous parametric downconversion (SPDC). Energy conserva- 
tion requires that the two downconverted photons have frequencies 
iO± = LOQ where 2a;o is the pump frequency. We again assume 
a Gaussian distribution of frequencies, centered around ujq, of the 
form e~ 2'''"±~'^''' = e~'2^ . We follow essentially the same cal- 
culational procedure as before, except now we enforce the require- 
ment that the incoming photon frequencies satisfy u)\+ u)2 = 2a;o . 
In this section we impose this condition in a manner that will allow 
us to obtain analytic solutions for the output probabilities. This 
will serve us as a simplified version of SPDC, and we will see in 
section [VIII that this model seems to give an upper bound to the 
mutual information obtained from a more realistic model of SPDC. 
The input state in this model is taken to be of the form 



dc.2e-t[(-i-»)'+("2-"»)'l |20>.= (/— / / d.e-" /(e)ai^aijO>, (44) 

V 7ry_oo J —00 



|20>^ |02><, = 

V 

the output probabilities are 

P(2,0|</.o) = P(0,2|</.o) 



1 1 - 

4 n 



2^ 



1 1 - 

- <^ 1 e 

2 I ri 



(41) 



(42) 



(43) 



In the absence of dispersion (a = /3 = 0) or in the narrow band- 
width limit {cr — > 00), we see that the coincidence rate P{l,l\{j>o) 
vanishes, while the other two probabilities are both equal to i . 

Note that there is no dependence on (f>o- We will see later that 
this fact manifests itself in a vanishing mutual information. 



where now ui-^- = luq + Q and ui— = uiq — Q + e. We can choose 
/(e) to be any function sharply peaked at zero with normalized 
integral (unit area under its graph). We then compute the output 
probabilities according to 



P{N^,Na\M = J dwiduJ2 |(A^c, Afd|Vm>P , (45) 



or equivalently, by applying the projection operators 



i,{N,,Nj) = J dull duj2\Nc,Na){N^,Na\. (46) 



B. Two-Photon Input with Anticorrelated 
Energies: Simplified SPDC model 

Case F: Simplified SPDC Foclc states. We now examine a 
case with two photons incident on the same input port and anticor- 
related in energy. We do this in the context of a simplified model of 



The auxiliary function f\{€) is necessary in this model in order 
to impose the constraint uii + lu2 = 2u>q without causing squares 
of delta functions to arise in the probability calculations. A more 
correct treatment of SPDC will follow in section [VIII 

The measurement outcomes, integrated over final frequency, are 
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then given by 
P(2,0|<Ao) = 



P(0,2|</.o) 



2 + e- 



1 / 6*1 

— cos [2<po + — 
ri V 2 



(47) 



4 

2 + e"~2S— 



cos 



6*2 a^^^S 



16t-20"2 



1/6*1 
— cos 2(/)o + — - 
n V 2 



(48) 



Vr2 



^ COS 



4>o ■ 



6*2 a^PL^ 



2 16r2o-2 



1 / Si 

— cos 2flio + — 
n V 2 



(49) 



As /3 increases, ri and r2 increase, leading to decreased visibility 
of all of the oscillating terms. 

Note also that in the case of zero dispersion (« = /3 = 0), the 
exact expressions for energy-uncorrelated (Case C, section 4.1) and 
energy-anticorrelated (downconverted) Fock states (Case F) are 
identical to each other. However the probabilities begin to diverge 
when dispersion is turned on. The same effect will be seen to occur 
for the uncorrelated and anticorrelated NOON states in the HOM 
interferometer (cases J and K, below). 

Case G: Simplified SPDC NOON states. Now the input 
state is taken to be a NOON state, -7= {|20)a,ct + |02>a,ct}. We find 
the measurement outcomes to be: 



P(2,0|</.o) = P(0,2) 



^'(1, 1100) 



4 

,2r2 



l + e" 



1 - e~ 



(50) 
(51) 



As in the uncorrelated case, the probabilities show no depen- 
dence on <j>0, and so have vanishing mutual information. In this case 
we also see that there is no dependence on the 2nd order dispersion 
coefficient /3. 

It is interesting to note what happens if we shift the phase of 
the photons in one input port by ^ before they hit the first beam- 
splitter. The input to the interferometer is now proportional to 
|2, 0} — |0, 2) . In this case, the interference in <j>0 reemerges, and 
the result is independent of a instead of /3. In fact, the counting 
probabilities turn out to be very similar to those of the NOON state 
incident on an HOM interferometer presented in the next section 
(case K). Moreover, these two cases have identical values for the 
mutual information. 

Case H: Simplified SPDC dual Fock state. The frequency- 
anticorrelated dual Fock input state 



gives the results 



de e 



P(2,0|</.o) = P(0,2|</.o) 



200- 



(52) 



(53) 



(54) 



V. DISPERSIVE HONG-OU-MANDEL 
INTERFEROMETER WITH NOON INPUT 

An alternative setup has been proposed to improve phase 
resolution(0|). In this section we examine this alternate version 
and compare it to the previous results. 

In this version, it is assumed that the NOON state is created in- 
side the interferometer, rather than at the input ports. Effectively, 
we need to remove the first beam splitter from the interferometer 



porta ^ pjt 



input 




output 



FIG. 3: Hong-Ou-Mandel interferometer with dispersion and 
non-dispersive phase shift (po in one arm. 



and use the NOON state as input to the remaining beamsplitter, 
which now acts as a Hong-Ou-Mandel (HOM) interferometer Il3ll . 
The setup is shown in figure |3] We again assume dispersion and 
phase-shift <j>o along one of the lines entering the beam splitter, 
neglecting absorption. Ignoring an overall constant phase of e*''"^, 
the scattering matrix now has the form 



1 / e»"*(") 



-1 



(55) 



■ ibu 



(56) 
(57) 



where 4>{ii>) is again given by equation 1121 Thus, 
, _ 1 

Coj — — 

V2 

dui = — ^ 
V2 

Case I: Single-photon NOON state in HOM interferome- 
ter. By the same methods as before, we can compute the counting 
rates for a given NOON state input. The one-photon NOON input 
state is 

(58) 



W=cl dc^e-f (ai+bi)|0>, 



for which we find 



Pii,o\M 




/ri 

^2 r 2 



4>o - 



^/rT 



2 4r2cr2 



(59) 



(60) 



Finally, assuming two-photons, with no correlation or with an- 
ticorrelation, we arrive at two additional cases (J and K). 

Case J: Energy-uncorrelated two-photon NOON state in 
HOM interferometer. 

The input state is 

(61) 

From this state, we arrive at the results: 
P(2,0|</.o) = P(0,2|flio) 



P(l,l|0o) 



1 ■ 



1-1- 



ri 

„2r2 



COS 2(t>o + 01 



cos 24)0 + 01 



2rfa^ 



2rfa2 



(62) 



(63) 
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a 

FIG. 4: (color online). Mutual information versus alpha for 
single photon cases (cases A, B, and I), plotted for the values 
f} = 0, a — 1. The mutual information is the same for all 
three cases, (a, /3, and a are in units of L~^uJq^, L'^^loq'^, 
and lJq^, respectively.) 




0.5 1 1.5 2 2.5 3 3.5 4 



P 

FIG. 5: (color online). Mutual information versus /3 for single 
photon cases (cases A, B, and I), for the values a — .5, a = 
1. (a, /9, and a are in units of L'^lJq^, L~^uJq^, and lJq'^, 
respectively.) 



Case K: Simplified SPDC two-photon NOON state in 
HOM interferometer. 

For the input 



\^) = C j 



do. e 



we compute: 

P(2,0|</.o) = P(0,2|</.o) 



^^(1,1100) 



1 / 6»i 

1 cos 200 H 

Vn" V 2 

1 / 6*1 

1 + —— cos [2<t>Q + — 
\/ri \ 2 



(64) 



(65) 



(66) 



In this last case, the results turn out to be independent of the 
first order dispersion coefficient, a. 



VI. COMPARISON AND DISCUSSION OF 
CASES A TO K 



0.25 
0.2 
0.15 



FIG. 6: (color online). Mutual information versus squared 
inverse bandwidth a for single photon cases (cases A, B, and 
I), for the values a = 1, /3 = .1. [a, (3, and a are in units of 
L^^ujq^, L~^ujq^, and uJq^, respectively.) 



The detection probabilities of the previous sections can now be 
combined with equation [2] to compute the mutual information for 
each of the experimental setups and inputs states. Plotting the 
results as functions of a, /3, and a, we find the results in figures 
l4l6l for single photon input and figures fTIQI for two photons. aL is 
given in units of i^q^ , while fSL and a are in units of i^q'^ ■ In the 
dispersionless limit, o, /3 — > 0, we find Shannon mutual information 
values that agree with those previously calculated in jij]. 

Only positive values of fS were graphed. However, the formulas 
of the previous sections work equally in the anomalous dispersion 
(negative /3) region. 

Note also that the four parameters a, /3, cr, L appear in all equa- 
tions only through the dimensionless quantities 



Ai 



/3L 



and 



A2 



(67) 



Thus, other parameter ranges can easily be obtained from those 
graphed here via appropriate rescaling of variables with the dimen- 
sionless ratios held fixed. 

A few conclusions are immediately clear from these graphs and 
from the equations of the previous sections, (i) First, the dual Fock 
states entering the Mach-Zehnder interferometer give identical re- 
sults as the NOON states entering the Hong-Ou Mandel interferom- 
eter (compare equations 1541 and 1541 to 1651 and 1661 or compare 1391 




FIG. 7: (color online). Mutual information versus alpha for 
two-photon cases, for the values a — 1, P — 0. Cases E and 
G vanish identically. Cases J and D are identical, as are cases 
H and K. (a, /3, and a are in units of L~^ujq^, L~^ujq'^, and 
loq'^, respectively.) 
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0.5 1 1.5 



3 3.5 4 4.5 5 



FIG. 8: (color online). Mutual information versus beta for 
two-photon cases, for the values a = 1, a — .5. (a, /3, and a 
are in units of L~^(jJq^, L~^uJq^, and ujq^, respectively.) 




FIG. 9: (color online). Mutual information versus squared 
inverse bandwidth a for two-photon cases, for the values a = 
1,13 = .1. (of, P, and g are in units of L^^uJq^ , L^^lJq^, and 
loq'^, respectively.) 



and 1401 to 1621 and I63I I. This is to be expected, since for N = 2 
the two cases are equivalent: the first beam splitter in the Mach- 
Zehnder interferometer turns a dual Fock input state into a NOON 
output state, which then strikes the second beamsplitter. The sec- 
ond beamsplitter can then be viewed as an HOM interferometer. 
Thus cases J and D are equivalent, as are cases H and K. (This 
equivalence will not for N > 2.) (ii) Second, the single-photon 
cases (cases A, B, and I) all give identical curves for the mutual 
information as functions of a, /3, and cr. The explanation for this 
is clear if the action of the first beam splitter on the input is exam- 
ined. Cases B and I are equivalent for the same reason mentioned 
in the previous point: they both lead to a one-particle NOON state 
in the portion of the interferometer before the dispersive element 
is reached, and so give the same output. Meanwhile, in case A, the 
output of the first beam splitter is the state (|01} -I- «|10)) ; this 
is similar to a NOON state, except one term is shifted in phase by ^ 
relative to the other. This converts the sines in the probabilities of 
cases B and I into the cosines of case A (equations 1251 and 1261 , but 
has no other efi'ect. Since the mutual information involves integrals 
from —IT to TT, interchanging sines and cosines inside the integrals 
has no effect on the mutual information. Unsurprisingly, the single 
photon cases generally result in lower mutual information than the 
two-photon cases, (iii) We see from the graphs for the two-photon 
states that the energy-uncorrelated and energy-anticorrelated ver- 
sion of each input give identical results for zero dispersion or zero 



bandwidth (cr = oo); however, the uncorrelated versions all drop ofi' 
rapidly to zero fidelity as the dispersion increases, whereas the anti- 
correlated (downconverted) input leads to a much slower drop, (iv) 
Two-photon NOON states incident on the MZ interferometer (cases 
E and G) have zero mutual information as anticipated earlier, (v) 
For fixed bandwidth and fixed quadratic dispersion coefficient /3, 
the two-photon downconverted NOON state in the HOM interferom- 
eter (case K) is independent of the linear coefficient a. However, 
it decays rapidly with increasing /3. (vi) Overall, the simplified 
SPDC-generated Fock states (case F) seem to hold up best in the 
presence of dispersion. This case starts off with a higher value of 
H at zero dispersion and decays more slowly as a and /3 increase. 
The only exception to this statement is when /3 is small, in which 
case the anticorrelated HOM NOON state (case K) works better at 
large a. 

A bit of insight into some of the properties of the 2-photon results 
may be obtained by considering the exponential decay factor 



(68) 



Ai and A2 are the dimensionless quantities defined in equation 1671 
In frequency-uncorrelated cases such as cases C and D, all of the 
00-dependent terms are multiplied by a factor of which arises 
from interference between e*'^("i' and 6**^^"^' terms, where uii and 
LJ2 are the frequencies of the photons entering the input ports. 
The relevant term is of the form e'[<*("i)+'^("2)l . As o ^ 00 or 
(7 — ► 0, we find that A2 — > and ^ 0, so that only constant 
(</>0-independent) terms survive in the limit. Thus, for large a 
or small cr, the dependence of the probability distributions on ipo 
decays exponentially, causing the mutual information to also decay 
rapidly. 

In contrast, for the frequency-anticorrelated cases, such as F and 
H, the term e'{"*("i)+'?'{"2)) becomes 



(69) 



with the ct-dependence cancelling. As a result, (/iQ-dependent terms 
occur without the exponentially decaying factor, allowing much 
slower decay of H at large a (or even no decay at all, as in case 
H). The slower decay at large dispersion is therefore a direct con- 
sequence of the quantum-mechanical correlations present in the 
initial state. 

As for the /3 dependence, we see that as /3 becomes large, both 

A2 

and ri become /3 independent, with C ~* £ * E'rid ri 1; thus 
all the curves approach constant values at large fS, with slopes 
of comparable order of magnitude. 

We turn now to one additional case, that of more realistic SPDC 
photon pairs, which we then proceed to compare with the simplified 
SPDC model already examined. 



VII. CASE L: SPDC 

Now we present results for the mutual information using a more 
realistic model for the parametric downconversion process. Numer- 
ically, the results turn out qualitatively (and for some parameter 
ranges quantitatively as well) to be very similar to those of the sim- 
plified SPDC model in the previous section; however we no longer 
will be able to present explicit analytic expressions for the mea- 
surement outcomes. 

There are many possible cases that could be considered, but we 
restrict ourselves here to the single case of coUinear type-II SPDC in 
a nonlinear crystal, with both of the outgoing photons entering port 
A of the dispersive Mach-Zehnder interferometer. We now have to 
consider the parameters of both the interferometer and the crystal. 
We allow the pump frequency to vary around central frequency 2u>o , 
with the deviation from the center of the distribution represented 
by 2Qp; in other words, the pump frequency is represented as 



ijjp = 2{ijJo + Up). 
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We once again assume a Gaussian distribution of frequencies, in this 

case represented by a weighting factor e~'2^'^p~^'^'^"^ = e~'^'^^p . 
The signal and idler frequencies are then 

uJs = ^ + n = a;o + f^p + n (70) 
uii = — f2 = a;o + — f^i (71) 

with ll)s + uji = Ljp. Suppose that the crystal is cut so that exact 
phase matching occurs at the central frequency 

fcp(2(^o) = fcs('^o) + fci(i^o)- (72) 

Then, assuming that terms quadratic and higher in the frequencies 
are small, the phase matching condition for the crystal gives us a 
condition on the wave-vectors of the form 14] 

Ak = kp(u)p) — ks(u)s) — ki{LUi) = ApSlp + ASl, (73) 



LU and a;', we have 

^=(2,0100) 



P(O,2|0o) 



where 4>{uj) is as defined in equation 1121 Note that 
$ ^ " ~" , ujQ — ^ depends only on the crystal properties, 

while 0(tt') depends only on the properties of the interferometer. 
The integrands of P(0, 2\{f>o) and P(2,0|</<o) factor in their depen- 
dence on these two sets of parameters; that of P{1, ll^o) does not, 
indicating the entangled nature of the |11) state. 

Given these output probabilities, the mutual information can 
once again be computed. In contrast to the previous sections, the 
analytic forms of the probabilities are too complicated to be enlight- 
ening, so we proceed to numerical calculations. Some examples are 
graphed in figures [T0l to ll4l The plots are expressed in terms of the 
new parameters b = ApLc and A = -4-. 

Examples of the dependence of H on the parameters of the pump 
beam (a), interferometer (a, /3), and nonlinear crystal (A, b) are 
given in figures 1101 through 1141 We see that, although H decays 
overall with increasing values of the dispersion paramaters in the 
interferometer, a and /3, there are oscillations superimposed on the 
decay, which are especially noticeable at low values of a and /3. 
This effect was in fact also present in the simplified SPDC model 
of the previous sections, but in the latter case the oscillations were 
too weak to be visible on the graphs. We see also that as either 
b or A increases (or equivalently, as A or Ap increases), the plots 
approach those of the simplified SPDG model. Since A and b are 
proportional to the crystal length, this means that the simplified 
SPDC model is an increasingly better approximation to real SPDC 
for longer crystals. It also appears from the numerical simulations 
that for a given set of parameter values a, /3, L, and a, the simpli- 
fied SPDC model provides an upper bound to H for the real SPDC 
cases with the same parameter values. The maximum information 
content clearly occurs for low dispersion in the interferometer, long 
nonlinear crystals, and large mismatch at uiq between signal and 



where Ap = 2k'p(2LJo) ~ k'^{u}o) — k[(uJo) and A = k'^(uJo) — k'g(ujQ). 

The wavefunction for the biphoton state entering the interfer- 
ometer is now 

|,^,„> = I dQdnp<f{np,n)al^^^^^^dl^^,^^_^\o), (74) 

where 

with normalization constant Af. Here, is the length of the non- 
linear crystal. Using this wavefunction, we can compute output 
probabilities as before. Denoting the frequencies at the detectors 

by 



(76) 



(77) 



(78) 



idler inverse group velocities in the crystal (large A). 



VIII. CONCLUSIONS 

In this paper, we have examined the effect of dispersion on the 
mutual information that interferometric photon-detection measure- 
ments carry about phase shifts. We have looked at a number of 
different situations involving two interferometer set-ups and sev- 
eral different types of non-classical input states. Comparing the 
results, we now have a precise and quantitative means to mea- 
sure the relative merits of different input states for various input- 
parameter ranges. As a by-product, we have shown that in some 
circumstances, parametric downconversion can be approximated by 
a much simpler model that is amenable to exact analytical analysis. 

Returning to the original question of which input state yields 
the most information about the phase shift, the graphs of the pre- 
vious sections yield fairly clear results. Restricting discussion to 
MZ interferometers for simplicity, we can see that for quantum 
interferometry in the presence of dispersion the entangled photon 
pair produced by downconversion has a clear advantage over other 
cases when input to a single port (Fock state input). This advan- 
tage does not exist in the case of an dispersionless interferometer, 
in which case the presence or absence of frequency correlations be- 
comes irrelevant for the information content. The only situation 
we have found in which another input is superior to the frequency- 
anticorrelated Fock input is when a is large but (3 small, in which 
case the anticorrelated dual Fock input is superior. These conclu- 
sions all hold when the simplified downconversion model of sec- 
tion |IV2] is a good approximation; the results of section FVIII implv 
that such conclusions weaken as the crystal becomes shorter. 



/ 



du dw' 
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FIG. 10: (color online). Mutual information versus squared 
inverse bandwidth a for SPDC. {a, {3, and o are in units of 
L~^uJq^, L~^u!q'^, and ujq^, respectively, b is in units oi ujq^, 
while A is dimensionless.) 



FIG. 13: (color online). Mutual information versus lambda 
(A — A/Ap) for SPDC. (a, and a are in units of L~^uJq^, 
L~^ujq^, and lJq'^, respectively, b is in units of lJq'^, while A 
is dimensionless.) 



Simplllied model (Case F] 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



FIG. 11: (color online). Mutual information versus alpha for 
SPDC with (J = 1, /3 = .1. (a, /3, and a axe in units of 
L~^liJq^ , L~^LLig'^ , and cOq^, respectively, b is in units of lOq^, 
while A is dimensionless.) 



Simplilied model (case F) 
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FIG. 12: (color online). Mutual information versus beta for 
SPDC with a = 1, a — .3. (a, /3, and a are in units of 
L~^ujq^, L~^uj^^, and ujq^, respectively. 6 is in units of luq'^, 
while A is dimensionless.) 
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FIG. 14: (color online). Mutual information versus b (6 



%Ap) for SPDC. 



(a, /3, and a are in units of L 







and luq , respectively, b is in units of luq , while 



A is dimensionless. 
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